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Abstract. The phenomenological analysis of the LHC data on px-spectrum and elliptic flow 
of J/'i!) and D mesons in PbPb collisions at ^snn = 2.76 TeV is presented. The charmed 
meson production pattern in PbPb collisions may be reproduced by two-component model 
HYDJET-b-l- including thermal and non-thermal components. The significant part of D-mesons 
is found to be in a kinetic equilibrium with the created medium, while J/tp-mesons are not. 


1. Introduction 

Studying the charmed hadron production is a particularly useful tool to probe properties of 
hot and dense matter created in relativistic heavy ion collisions. The modern pattern of multi¬ 
particle production in (most central) heavy ion collisions at RHIC and LHC agrees with the 
formation of hot strongly-interacting matter with hydrodynamical properties (“quark-gluon 
fluid”), which absorbs energetic quarks and gluons due to their multiple scattering and medium- 
induced energy loss (see, e.g., m)- Within such paradigm, a number of questions on heavy 
flavor production arises. Are heavy quarks thermalized in quark-gluon plasma? What is the 
mass dependence of in-medium quark energy loss? Are charmed hadrons in a kinetic equilibrium 
with the medium? How the specific pattern of quarkonium suppression is related to the interplay 
between thermal and non-thermal mechanisms of J/^-meson production? In this paper, some 
results of the phenomenological analysis of LHC data on momentum spectra and elliptic flow 
of charmed hadrons (D, J/V^) in PbPb collisions at y^SNN = 2.76 TeV within two-component 
model HYDJET-|--|- are presented. The comparison with RHIC results is also discussed. 

2. HYDJETH—[- model 

HYDJET-|--|- is the model of relativistic heavy ion collisions, which includes two independent 
components: the soft hydro-type state (“thermal” component) and the hard state resulting from 
the medium-modified multi-parton fragmentation (“non-thermal” component) [2]. 

The soft component represents the hadronic state generated on the chemical and thermal 
freeze-out hypersurfaces obtained from the parametrization of relativistic hydrodynamics with 
preset freeze-out conditions (the adapted event generator FAST MC [Sill]). Hadron multiplicities 



are calculated using the effective thermal volume approximation and Poisson multiplicity 
distribution around its mean value, which is supposed to be proportional to a number of 
participating nucleons for a given impact parameter of an AA collision. To simulate the elliptic 
flow effect, the hydro-inspired parametrization is implemented for the momentum and spatial 
anisotropy of a thermal hadron emission source. 

The approach used for the hard component is based on the PYQUEN partonic energy 
loss model [S]. The simulation of single hard NN sub-collision by PYQUEN is constructed 
as a modification of the jet event obtained with PYTHIA_6.4 generator [B]. It takes into 
account medium-induced rescattering, collisional and radiative energy loss of hard partons in 
the expanding quark-gluon fluid, realistic nuclear geometry and nuclear shadowing effect. The 
mean number of jets produced in an AA event is calculated as a product of the number of binary 
NN sub-collisions at a given impact parameter per the integral cross section of the hard process 
with the minimum transverse momentum transfer The partons produced in hard processes 
with the momentum transfer lower than p™", are considered as being “thermalized”. So, their 
hadronization products are included “automatically” in the soft component of the event. 

The input parameters of the model for soft and hard components have been tuned from fitting 
to heavy ion data on various observables for inclusive hadrons at RHIC [2] and LHC [?]■ 

Charmed meson production in HYDJET-|--|- includes soft and hard components as well. 
Thermal production of D, J/'ijj and Ac hadrons is treated within the statistical hadronization 
approach B Momentum spectra of charm hadrons are computed according to the 

thermal distribution, and the multiplicities Nc {C = D, J/'ll), Ac) are calculated through the 
corresponding thermal numbers as Nc = 7c‘'A"c^, where 7 c is the charm enhancement factor 
(or charm fugacity), and Uc is the number of charm quarks in a hadron C. The fugacity 7c can 
be treated as a free parameter of the model, or calculated through the number of charm quark 
pairs obtained from PYTHIA and multiplied by the number of NN sub-collisions. Non-thermal 
charmed hadrons are generated by PYQUEN taking into account in-medium energy loss of heavy 
(b, c) quarks within “dead-cone” generalization [10] of BDMPS model [111 112] for radiative loss 
and high-momentum transfer limit daEiiis] for collisional loss. 

3. J/V'-meson production in PbPb collisions at the LHC 

Some time ago it was demonstrated in [16] that if J/V’-mesons are produced in HYDJET-I—|- 
at the same freeze-out parameters as for inclusive (light) hadrons, then simulated px- and y- 
spectra in 20% of most central AuAu collisions at RHIC energy = 200 GeV are much wider 

than ones measured by PHENIX m- However the assumption that the thermal freeze-out for 
J/'0-mesons happens at the same temperature as chemical freeze-out (with reduced collective 
velocities) allows HYDJET-I—|- to reproduce the data. It turned out that the similar situation 
holds at the LHC. Figured] shows the comparison of HYDJET-I—|- simulations with the ALICE 
data [18] for px-spectrum of JZ-^-mesons in 20% of most central PbPb collisions at y^snn = 2.76 
TeV. One can see that if thermal freeze-out for J/ip^s happens at the same temperature as 
chemical freeze-out (with reduced collective velocities), then simulated spectra match the data 
up to px ~ 3 GeV/c. The discrepancy at high px may indicate on the necessity to tune PYTHIA 
for charmonium production in this kinematic domain. The px-dependence of the elliptic flow 
coefficient V 2 measured by ALICE [19] is also reproduced by HYDJET-I—|- simulations in such 
a case (figure [2|). The contribution of non-thermal component is important at high px- 

Thus we conclude that J/'ip-m.osoxi production in heavy ion collisions at RHIC and LHC 
may be reproduced by HYDJET-|--|- model under the assumption that thermal freeze-out of 
charmonia happens appreciably before thermal freeze-out of light hadrons [20], presumably 
at chemical freeze-out (with reduced radial and longitudinal collective velocities). Our 
interpretation of this is that the significant part of J/^Z-mesons (up to px ~ 3 GeV) is 
thermalized, but without being in a kinetic equilibrium with the medium. 
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Figure 1. Transverse momentum spectrum 
of J/?/;-mesons at rapidity 2.5 < 7/ < 4 in 20% 
of most central PbPb collisions at y^snn = 
2.76 TeV. The points are ALICE data [18], 
histograms are simulated HYDJET++ events 
(red - freeze-out parameters as for inclusive 
hadrons, blue - “early” thermal freeze-out). 


Figure 2. Elliptic flow coefficient V 2 {pt) of 
J/t/;- mesons at rapidity 2.5 < y < 4 for 20- 
40% centrality of PbPb collisions at = 

2.76 TeV. The points are ALICE data |19j . 
histograms are simulated HYDJET-|--|- events 
(blue - soft component, red - hard compo¬ 
nent, black - both components). 




Figure 3. Transverse momentum spectra of 
D-mesons at rapidity | y |< 0.5 in 20% of most 
central PbPb collisions at y'SNN = 2.76 TeV. 
The points are ALICE data [22|, histograms 
are simulated HYDJET-I—|- events (blue - soft 
component, red - hard component, black - 
both components). 


Figure 4. Elliptic flow coefficient V 2 {pt) of 
D'^-mesons at rapidity | y |< 0.8 for 10-30% 
centrality of PbPb collisions at \/snn ~ 

TeV. The points are ALICE data [23], his¬ 
tograms are simulated HYDJET-I--|- events 
(blue - soft component, red - hard compo¬ 
nent, black - both components). 


4. D-meson production in PbPb collisions at the LHC 

At first we have checked that px-spectrum of D-mesons measured by STAR [21] in 10% of most 
central AuAu collisions at y^SNN = 200 GeV is reproduced by HYDJET-I--|- simulation with 
fugacity 7 c = 7 and the same freeze-out parameters as for J/ 7 /^-mesons (but not for inclusive 
hadrons). It means that D-mesons like J/ 7 /^-mesons are not in a kinetic equilibrium with the 
medium at RHIC. However the situation gets changed dramatically at the LHC. Figure [3| shows 
the comparison of HYDJET-I—|- simulations with the ALICE data [22j for px-spectra of D*^, D^ 
and D*^ mesons in 20% of most central PbPb collisions at y^SNN = 2.76 TeV. The freeze-out 
parameters as for inclusive hadrons have been used, the fugacity value 7 c = 11.5 being fixed 
from J/?/) yield. The model results are close to the data within the experimental uncertainties. 
Note that px-spectra and nuclear modification factors 72 aa(pt) (the ratio of particle yields in 
AA and pp collisions normalized on the mean number of NN sub-collisions) of D-mesons are 
























































reproduced by HYDJET++ up to very high px- So the treatment of heavy quark energy loss 
in hard component of the model (PYQUEN) seems quite successful. In addition, HYDJET++ 
is able to reproduce the ALICE data [23] on px-dependence of the elliptic flow coefficient V 2 
(figure m) . 

Thus in contrast to RHIC, thermal freeze-out of D-mesons at the LHC happens 
simultaneously with thermal freeze-out of light hadrons. Therefore the signihcant part of D- 
mesons (up to px ~ 4 GeV/c) seems to be in a kinetic equilibrium with the medium. 

5. Summary 

Momentum spectra and elliptic flow of D and J/'ijj mesons in PbPb collisions at the LHC are 
reproduced by two-component model HYDJET-|--|- including thermal (soft) and non-thermal 
(hard) components. Thermal freeze-out of D-mesons happens simultaneously with thermal 
freeze-out of light hadrons. Thermal freeze-out of J/'0-mesons happens appreciably before freeze- 
out of light hadrons, presumably at chemical freeze-out (with reduced radial and longitudinal 
collective velocities). Thus the significant part of D-mesons (up to px ~ 4 GeV/c) seems to 
be in a kinetic equilibrium with the medium, while mesons are not. The production of 
charmed mesons at high transverse momenta is determined by non-thermal charm production 
mechanism and medium-induced heavy quark energy loss. 
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